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Abstract 

We calculate the Raman spectrum of the two-dimensional (2D) spin-1/2 
Heisenberg antiferromagnet by exact diagonalization and quantum Monte 
Carlo techniques on clusters of up to 144 sites and, on a 16-site cluster, by con- 
sidering the phonon-magnon interaction which leads to random fluctuations 
of the exchange integral. Results are in good agreement with experiments on 
various high-T c precursors, such as La2Cu04 and YBa2Cu30g.2- In particu- 
lar, our calculations reproduce the broad lineshape of the two-magnon peak, 
the asymmetry about its maximum, the existence of spectral weight at high 

energies, and the observation of nominally forbidden Ai g scattering. 
PACS numbers: 78.30.Hv, 75.50.Ee, 74.72.Jt 
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Introduction. — Raman scattering is a powerful technique to study electronic excitations 
in strongly correlated systems. Recently, much attention has been given to the anomalous 
magnetic scattering with a very broad and asymmetric lineshape observed in the Raman 
spectra of the parent insulating compounds of high-T c superconductors, such as La2Cu04, 
and YBa2Cu 3 6 . 2 at around 3230cm -1 and 3080cm -1 , respectively ]!]]. The selection rules 
associated with this peak are also anomalous. While the spin-pair excitations scatter pre- 
dominantly in the allowed B\ g channel, there is also a significant contribution in the nomi- 
nally forbidden A\ g configuration, as well as much weaker B2 9 and A2 9 scattering ]]]. 

Previous theoretical studies on the spin- 1/2 Heisenberg model for 2D square lattices have 
computed the Raman spectra and its moments for a nearest-neighbor interaction and 
only the moments when spin interactions along the plaquette diagonal were also included. || 
These show good agreement with experiments in regard to the position of the two-magnon 
peak, but they fail to account for the spectral shape, and its enhanced width. 

Several schemes have been considered to resolve this problem. Initially, from the analysis 
of the moments it was proposed that strong quantum fluctuations were responsible for the 
broadening (see, e.g., Ref. [0,§J). However, recent studies of spin-pair excitations in a spin- 
1 insulator, MPS3, show a width comparable to that of the spin-1/2 cuprates |7[]. This 
questions the view that the observed anomaly is due to large quantum fluctuations intrinsic 
to spin-1/2 systems. We remark that the measured widths are 3-4 times larger || than 
those predicted by Canali and Girvin || within spin-wave theory using the Dyson- Maleev 
transformation, even when processes involving up to four magnons are taken into account. 
The work by Canali and Girvin |5j and other groups P,|10f present convincing evidence 
that the observed anomalous features of the magnetic scattering cannot be satisfactorily 
explained by only considering quantum fluctuations. 

In order to explain the observed anomalously broad and asymmetric lineshapes, it seems 
then necessary to invoke an additional process. Here, we consider the interaction between 



magnon pairs and phonons |TI[ . This mechanism is motivated in part by recent experimental 



observations of a strong broadening of the B\ g and an enhancement of the A\ g scattering 



with increasing temperature |12[ . In our approach we consider the phonons as static lattice 
distortions which induce changes, SJij, in the exchange integral J of the undistorted lattice. 
We calculate the Raman spectra for a nearest-neighbor Heisenberg model using a nearest- 
neighbor Raman operator in the quenched-phonon approximation which, like the Born- 
Oppenheimer approach, focuses on the fast (high-energy) magnon modes and freezes the 
slow (low-energy) phonons. This approximation is valid for the cuprates because there is 
a clear separation of energies between the magnetic and vibrational modes. For instance, 
in YBa2Cus06 the characteristic Debye frequency is about 340cm _1 while the two-magnon 
excitation is ~ 3080cm -1 . 

Raman Lineshape without phonon-magnon coupling. — The isotropic Heisenberg Hamil- 
tonian is given by H = JJ2<ij> Sj • Sj , where the notation is standard, and only nearest 
neighbor interaction is assumed. For the cuprates, the exchange integral is J ~ 1450-K" ~ 
0.12eV. In our study, we obtained the ground state |0o) of H on finite 2D square clusters 
with iV spins and periodic boundary conditions using a Lanczos (N = 16, 26), and Quantum 
Monte Carlo (QMC) (N = 144) algorithms. We studied zero and finite temperature spectra 
associated with the nearest-neighbor scattering operator [1-4] 

R— ^ (Ei rac ■ Ojj)(E sc ■ Q~ij)Sj ■ Sj, (1) 

<ij> 

where E inCtSC corresponds to the electric field of the incident and scattered photons, and 5V,- 
is the unit vector connecting sites i and j. In the cuprates, and for nearest-neighbors only, 
the irreducible representations of R are B\ g , Ai g , and E. We concentrate mainly on the 
dominant B lg scattering, e.g., E inc oc x + y and E sc oc x — y. 
The spectrum of the scattering operator can be written as 

I(u) = £ \(<p n \R\<p Q )\ 2 5(u - (E n - E )), (2) 

n 

where <p n denotes the eigenvectors of the Heisenberg model with energy E n . When doing 
exact diagonalizations on small clusters, the dynamical spectrum is extracted from a 
continued fraction expansion of the quantity 



I(u) = --Im{(j> Q \R | n 1 ■ W R\<Pv), (3) 

where e is a small real number introduced in the calculation to shift the poles of Eq. 3 into 
the complex plane. In the QMC simulations, the imaginary-time correlator (R(t)R(O)) is 
calculated and I(oo) is obtained by numerically continuing this function to real frequencies 
using a maximum entropy procedure ]13| . 

Our calculated B\ g spectra are shown in Fig. 1(a). They were obtained from exact di- 
agonalization (N = 16) and QMC (N = 144) studies of the Heisenberg Hamiltonian on 
square lattices. The two-magnon excitation observed experimentally lies around 3J, which 
is in good agreement with the location of the main peak obtained from exact diagonaliza- 
tion in Fig. 1. The position of this peak can be understood in terms of the Ising model, 
which corresponds to the limit of the anisotropic Heisenberg Hamiltonian when no quantum 
fluctuations are present. In its ground state, the Ising spins align antiferromagnetically for 
J > 0. Within this model and for a 2D square lattice, the incoming light creates a local 
spin-pair flip at an energy 3J higher than the ground state energy. This argument remains 
approximately valid even in the presence of quantum fluctuations PH^]. Our results indi- 
cate that the two-magnon excitation is at 2.9757J, 3.0370J, and 3.2J for the 16-, 26-, and 
144-site square lattices, respectively. Finite-size effects are small because of the local nature 
of the Raman operator. For the 144-site lattice, the QMC calculation was carried out at a 
temperature T = J/4. The slight shift of the peak position, compared to the T = results 
for the smaller clusters, is consistent with the finite-T exact diagonalization results of Ref. 
0. Statistical errors, absent in the exact diagonalization results but unavoidable in any 
stochastic simulation, enhance the width of the 144-site spectrum. These results confirm 
that neither finite-size effects nor finite temperature can account for the discrepancies with 
the experimental spectra. 

Lineshape Anomaly. — The Raman spectra obtained from the pure Heisenberg model 
(see Fig. 1) shows good agreement with experiments in regard to the two-magnon peak 
position, but the calculated width is too small. We will consider here the coupling between 
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the magnon pair and phonons ||TT| , p!2| to account for the observed wide and asymmetric 
lineshape. Our mechanism relates to that proposed by Halley [|14]] to account for two-magnon 
infrared absorption in, e.g., M11F2. 

Quantum and thermal fluctuations distort the lattice. The exchange coupling, which 
depends on the instantaneous positions of the ions, can be expanded in terms of the their 
displacements from equilibrium u. Keeping only the dominant linear terms: Jij{r) = = 
J+5Jij = J+u- VJjj(R). Here, 5Jij represents the instantaneous value of u- VJjj(R), where 
R denotes the equilibrium position of the ion carrying the spin (located at r = R + u) . In 
the quenched-disorder approximation, the effective Hamiltonian is 

Hi= ^(J+SJ^-Sj, (4) 

<ij> 

where \SJij\ < J is a random variable corresponding to taking a snapshot of the lattice. 
This new Hamiltonian is no longer translational invariant. 

In our study, the random couplings 5Jij were drawn from a Gaussian distribution 
P(5Jij) = exp (—(5 Jij) 2 /2a)/ \p2m a. I(uj) was obtained as the quenched average over 
m ~ 1000 realizations of the randomly distorted lattice. The quenched average of an oper- 
ator O is defined by ((O)) = ^ X^i(0o(j)|O|0o(j)), where <f>o(j) is the ground state of the 
jth realization of the disordered system. 

In Fig. 1(b) we show the B\ g Raman spectrum from Eq. (1) for a 16-sites square lattice 
with a ~ 0.4 J, which we found to agree best with experimental spectra [0. Our calculations 
do not consider the effect of frozen phonons on the scattering operator R. Notice that the 
coefficients pertaining to R are generally unrelated to the matrix elements of the system's 
Hamiltonian (e.g., dJ/dQ in H bears on e 2 /r, while the corresponding terms oc QSiSj in R 
bear on the dipole moment). In particular, and unlike the case without phonons, the fully 
symmetric Ai g component of the scattering operator does not commute with H. 

We find that the three main features observed in the B\ g configuration |IJ, namely, 
the broad lineshape of the two-magnon peak, the asymmetry about its maximum, and the 
existence of spectral weight up to w ~ 7 J are well reproduced. Beyond the two-magnon 



peak, there is a continuum of phonon-multi-magnon excitations. The small feature around 
uj ~ 5.5J (for < a < 0.3J) is compatible with a four-magnon excitation. 

Magnetostriction. — Since the effects of the phonon-magnon interaction (i.e., magne- 
tostriction) have not been extensively studied by theoretical work in the cuprates, a few 
comments are in order. The coupling between the spin and strain degrees of freedom 
modifies both elastic and magnetic properties. In fact, there are extensive studies on the 
(sometimes very strong) influence of elasticity on magnetism []15|-|17|. Mattis and Shultz |15 



considered the influence of uniform compression (i.e., all bonds equally distorted) in their 
classic study of magnetothermomechanics. Their results were criticized |I| for ignoring the 
effects of phonons (i.e., local fluctuations in the bond lengths, which are taken into account 
in the present work). Recently, giant magnetostrictive effects have been reported in several 
high-T c superconductors ||18|| . Also, important magnetostrictive effects have been reported 



in heavy-fermion ]19[ and low-T c [^0] superconductors. 



Superexchange-Phonon Coupling. — The width of the Gaussian distribution, a, represents 
changes in J due to large incoherent atomic displacements. Thus, one can write o ~ 
\(5\nJ/5Q)(Q)\ where (Q) is an average zero point motion (at T = 0) and (5J/5Q) is a 
weighted average of V Jy with respect to the displacement of all the ions participating in the 
exchange. Parenthetically, it is trivial to treat the case T ^ by increasing 5 J. Let r be the 
Cu-Cu distance, v the sound velocity, and M an effective reduced mass for the ions. A simple 
calculation gives (Q) jr ~ [Mvr /h)^ 1 ^ 2 ~ 0.05 which is consistent with X-ray measurements 
of the mean displacement of oxigen atoms normal to the layers [pi ,p2 ] . While VJj-, is not 



known for most phonons, values for longitudinal acoustic modes can be gained from the 



r-dependence of J in the form J(r) ~ r~ a or <91nJ/<91nr = —a [23]. For conventional 
transition metal oxides and halides, 10 < a < 14 p3"] , in reasonable agreement with the 
theoretical estimate a = 14 . For the cuprates, high-pressure Raman measurements [25 



and material trends [[2(J give, respectively, a ~ 5 — 7 and a ~ 2 — 6. These values translate 
into cr « (0.1 — 0.35) J. We emphasize that the relevant incoherent <5J's (or 5Q J s) of our 
case are much larger than those in pressure studies involving coherent motion of ions (see, 



e.g., the discussion in p. 466 of fll4l). Thus, we must use larger a (a ~ 0.4 J). 

Finally, we would like to stress that not every kind of disorder gives rise to the observed 
broadening of the spectrum. For instance, disorder by point defects or twinning planes will 
not produce such an effect. Also, it is observed in experiments that the Raman linewidth 
broadens with increasing temperature [T^|. This is a strong indication of a phonon mecha- 



nism for the broadening. 

A\ g and B2 g Symmetries. — For the A\ g symmetry, the undistorted Raman operator 
commutes with the Heisenberg Hamiltonian, and no scattering can take place. However, 
the addition of disorder changes the commutator and can produce an A\ g signal. Instead, the 
silent B 2g channel remains forbidden within our nearest-neighbor Raman operator. Fig. 1(c) 
shows the comparison between our numerically obtained A lg spectra (for a ~ 0.4J) and the 
experimental results The agreement between theory and experiments is reasonably 

good. We stress that the A\ g scattering follows naturally from our model unlike approaches 
relying on additional hypotheses, like, for instance, diagonal-nearest-neighbor couplings ||, 
4-spin terms [[J, new fermionic quasiparticles [|27| , or spinons. For a detailed discussion of 
these and other proposed explanations of the lineshape anomaly, see . 

Extensions. — The mid-infrared optical absorption in undoped lamellar copper oxides 
show broad features which are believed to originate from exciton-magnon absorption pro- 
cesses 



2S . Instead, these results could be interpreted as due to the interaction between 



phonons and magnons [p9| |. Also, excellent agreement with conductivity experiments has 



been recently achieved by the inclusion of phonon- induced strong-disorder |p0[ . 

Summary. — We find that light scattering spectra by spin excitations is caused by intrin- 
sic spin-spin interactions and by interactions with phonons. We provide strong evidence that 
the two-magnon Raman peak is strongly modified by coupling to low-energy phonons which 
randomly distort the lattice. Our calculations are in good agreement with experiments and 
provide a simple explanation of four puzzling features of the data: the broad lineshape of 
the two magnon peak, the asymmetry about its maximum, the existence of a spectral weight 
at high energies, and the observation of nominally forbidden A\ g scattering. 
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FIGURES 

FIG. 1. Normalized Raman cross section, I (uj) / I max , versus u/J, for the spin-1/2 Heisenberg 
model with iV sites, (a) B\ g Raman spectra obtained from exact diagonalization with N = 16, 
T = 0, and e = 0.1J; and from QMC with N = 144 and T = J/4. B lg (b) and A lg (c) 
spectra obtained from exact diagonalization (N = 16) with randomness in the exchange integral 
representing the interaction between spin-pairs and the phonons. The continuous, dashed and 
dotted lines in (b) and (c) correspond, respectively, to a = 0.3J, 0.4J, and 0.5 J. For comparison, 
the experimental results (from Ref.[6]) are shown. 
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